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Abstract—Comprehensive Assembly and Test Thermo 
Vacuum Chamber (CATVAC) facility is used _ for 
performing thermo vacuum performance test and thermal 
cycling tests on spacecraft and spacecraft appendages. The 
thermal system consists of a thermal shroud as an active 
element and test object as the passive one. Thermal cycling 
test profiles are temperature profiles of ramp and dwell on 
the test object with the process value control channel on it. 
The failure of the control system input thermocouple on the 
test object can disturb the control process. The designed soft 
sensor can infer the temperature of the control thermocouple 
channel from other measurements and provide redundancy 
to the physical channel. Archived thermal data was used for 
generating a process model with measurement channels as 
input and control channel value as output. The model input 
data was selected by studying their correlation to existing 
control channel measurements. The selected data was 
checked for validity using statistical method. The model 
output was obtained with another set of input data and 
validated by comparison with actual sensor values. The 
result shows excellent matching of soft sensor values to that 
of actual sensor. 


Keywords—ARX Model, PCA, Soft Sensors, System 
Identification, Thermal cycling. 


I. INTRODUCTION 


Comprehensive Assembly and Test Thermo-Vacuum 
Chamber (CATVAC) thermo-vacuum facility is used for 
testing space-crafts and their appendages in the space 
simulated thermal and vacuum environment. The facility 
is capable of simulation of vacuum of the order of 10° 
mbar and temperatures from -180°C to +120°C. 
Temperature control inside the chamber is met using 
liquid or gaseous nitrogen as media using a closed loop 
control system. The thermal control system controls the 
flow as well as temperature of the nitrogen based on set- 
point and the process value of the control channel on the 
test object. Thermal Cycling Facility (TCF) an additional 
attachment to the CATVAC chamber is used for rapid 
thermal cycling of appendages after being placed inside 
the chamber. The appendages could be dual gridded 
reflectors, antennae and solar panels of spacecrafts. Since 
the solar panels of the spacecrafts are exposed to the cold 
Space environment without any protection, — their 
qualification to those temperatures before integration to 
the satellite is must. The paper describes the application 
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Plant data selection from archives 
Model input selection and outlier 
detection 
Model structure selection 


Model Identification 
Model Validation 





Implementation 


Fig. 1. Schematic for design of soft sensor 


of soft sensor to the data obtained from the test of a 
dummy solar panel used for the qualification of the 
facility. The temperature data at different locations on the 
spacecraft is monitored using thermocouples connected to 
the data acquisition system. The control thermocouple 
which gives the temperature feedback to the control 
system is also placed on the solar panel. The PID control 
system accordingly varies the nitrogen density and 
temperature. In case of failure of the control channel or 
data blackout the control system will fluctuate or 
malfunction. Soft sensors or inferential sensors are an 
important tool to work comparatively with the physical 
sensors [1]-[5] so that even in case of any abnormal event 
of the physical sensor the reliability of the control system 
is not compromised. Historic data from the plant was used 
for designing the soft sensor. The available data from the 
process was analyzed using statistical methods and 
candidate inputs which have the maximum impact on the 
output temperature were selected. The selected data was 
filtered using techniques for outlier detection. A suitable 
model structure was selected and was validated with 
another set of data. The scheme used for the design of the 
soft sensor is shown in figure 1. 
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Fig. 2. Process flow and instrumentation diagram at CATVAC 


Ul. SETUP AND CONFIGURATION 


Solar panel is mounted inside TCF facing the shroud and 
thermally isolated from the mounting fixture as shown in 
figure 2. The TCF is then integrated into the chamber 
along with connections of the thermal line for flow of 
gaseous nitrogen inside the shrouds of the TCF. The 
chamber is then pulled to hard vacuum of 10° mbar range 
and start of thermal cycling subsequently. 20 numbers of 
thermocouples were placed on the solar panel facing the 
shroud and monitored through the data acquisition system. 
The control channel for which the soft sensor is designed 
is considered as the output of soft sensor. 


If]. INPUT SELECTION AND OUTLIER DETECTION 


The data acquisition system of the test setup acquires data 
from several thermocouples and stores it. The bubbled 
shroud through which nitrogen circulates is designed to 
maintain temperature uniformity throughout its surface. 
The data obtained from all the channels is shown in figure 
3. 


A. Input Selection 


In order to reduce the computational load for calculations 
the number of inputs was reduced to 6. Channels located 
physically away from the control channel were discarded 
from calculation. A simple correlation method [1] used to 
identify dependence of the control channel output on the 
inputs will be used. The scatter plot of each selected input 
against the output is shown in figure 4. Definite figure 
comes out of first four plots indicating strong linear 
correlation between the inputs and the outputs. Since plots 
with input 5 and input 6 did not indicate any definite 
relationship, they were discarded from __ further 
identification calculations without significant difference. 
Samples from the selected inputs were divided into two, 
Where around half samples were taken for model 
identification and other half were selected for validation 
of the identified model. Figure 5 shows the plot of the 
selected inputs 
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Fig. 3. Temperature data from all the channels 





Fig. 4. Scatter plot of selected input against output 
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Fig. 5. Plots of selected four channels 


B. Outlier Detection 


Outlier detection involves careful analysis of the obtained 
data and filtering of same for missing data, disturbances, 
and glitches due to malfunctioning of data acquisition or 


communication systems. Any set of data conveying 
improper information will have a negative impact on the 
quality of the process model. The input data obtained was 
first normalized for scaling the data. Statistical method for 
outlier detection reviewed in [9] was used. In this method 
Principal Component Analysis [10] is applied to the input 
variables and outliers are detected by plotting Jolliffe 
parameters and a sample is identified as outlier if its value 
exceeds the rest [1]. 
The co-variance matrix C can be computed using (1) from 
the input data in matrix X. 
Ca XX (1) 

The Principal Component decomposition of matrix X 
gives the matrix T given in (2). 

Xie (2) 
P contains the eigenvectors of the matrix C. The variables 
in T are called as scores. Jolliffee parameters [1] are 
computed using (3). 
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where p is the number of inputs, q is the number of 
principal components whose variance is less than one, Six 
is the i” sample of the k” principal component, 1, is the 
variance of the k" component. Plots of computed Jolliffe 
parameters are shown in figure 6. The inspection of 
Jolliffee parameters did not reveal any significant outliers 
in the data hence there was no need of extrapolation of 
missing or invalid data. 


IV. MODEL IDENTIFICATION 


Model of the system is the relation between the inputs and 
outputs of the system. Data driven approach of model 
identification to perform black-box modelling is used in 
this application [6]. The identified variables are given as 
input to the model and effect of these variables on the 
dependent variable or the output is studied for relation. 
The dynamics of the process indicates the suitability of 
parametric model structure to it; hence Autoregressive 
Models with Exogenous Inputs (ARX) was used [7]-[8]. 
The structure of the ARX model is shown in figure 7. 
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Fig. 6. Jolliffe Parameters (a) d,? (b) dy’ (c) d3” 





Fig. 7. Structure of ARX model 


where u(n), e(t) and y(n) represent the input, noise and 
output respectively. The noise error is neglected in 
calculations. The ARX model structure is given by (4). 


y(t)+ a, v(t = 1)+ eat Ang V(t a na) 
= hu(t—nk)+...+b,,u(t—nk —nb +1) 
where na is the number of poles, nb is the number of zeros 


and nk is dead time of the system. Another representation 
of the ARX model is given by (5)-(6). 


A(q)v(t)= B(q)u(t) + e(¢) (5) 
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where q'| is the delay operator. Using system 


identification techniques the model was identified. The 
method of least squares was used for identification of the 
ARX model. The identified parameters of the model are 
given by (7). 


A(q)=1-1.666q' +0.41487? +0.43397° 
~0.181 44 

Bq) =-0.4644q" +0.40797 * +0.40097° 
~0.345 44 

B,(q) = 0.72367 ' —0.67847? —0.520ly~* 
+0.4765q* 

B,(q) =—-0.16847! +0.15367° +0.287 4? 
~0.271974 

B,(q) =-0.016147! —0.00497Y? +0.11137° 
0.09024 


(7) 


V. MODEL VALIDATION AND RESULTS 


Once the model of the process is identified it has to be 
proved on other set of data for its robustness. A set of data 
from the measured samples was identified to be used for 
model validation to identify loss of model quality due to 
overfitting. This set of data was not used for identification 
of the model. The plot of the input data to be used for 
validation is shown in figure 8. This data was given as 
input to the model identified earlier and the output 
generated with the model was compared with the original 
output which was obtained from the physical channel. The 
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ig. 8. Plot of input data used for validation 


mydata; measured 
arx4412: fit: 99.81% 





Temperature (simulated and original) PC 


1000 2Q00, ples 2000 4000 


es 


ig. 9. Plots showing comparison of model output and physical channel 
data 
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Fig. 10. Difference of Soft Sensor and Actual Data 


result is shown in figure 9 where a perfect fit between 
model arx4412 output and actual obtained data is seen. 
Figure 10 shows the plot of the difference between the 
soft sensor simulated process values and actual sensor 
values. Thus the model proved its suitability for the soft 
sensor and can be implemented online as backup for 
physical channel. 


VI. CONCLUSION 


The designed soft sensor can be implemented using 
software platforms for virtual instrumentation [11]-[13]. 
These can be implemented using dynamic system 


identification where real time soft sensing can be done. 
The virtual instrumentation software can be used for real 
time data acquisition, pre-processing and __ for 
programming of soft sensor. The online system 
identification methods can capture entire range of system 
dynamics as the process values goes through different 
phases. The designed soft sensor will be used in 
maintaining reliability and robustness of the control 
system. The work will be further extended in design of 
soft sensing of nonlinear dynamics involved in 
temperature variations inside complex test objects like 
dual-gridded reflectors, antennae and other payloads 
subjected to thermal cycling. 
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